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Abstract 


In  this  report,  a  fast  algorithm  for  computing  scattered  fields  with  the 
use  of  physical  optics  (PO)  equivalent  approximation  in  half-space  is 
presented.  The  theoretical  basis  for  the  algorithm  and  the  derivation 
of  formulas  used  in  the  algorithm  is  presented.  The  algorithm  is  used 
to  compute  the  radar  cross  sections  (RCSs)  of  several  objects.  The  RCS 
of  the  objects  computed  by  the  algorithm  is  compared  with  those  that 
were  computed  with  the  method  of  moments  (MOM).  The  results 
presented  are  found  to  be  accurate  when  the  target  dimensions  are 
greater  than  or  equal  to  IX,  where  X  denotes  the  wavelength.  It  is 
concluded  that  the  PO  algorithm  presented  in  this  report  can  be  used 
for  majority  of  applications  as  it  captured  all  the  salient  features  of  the 
targets. 
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1.  Introduction 


In  general,  scattered  fields  from  complex  objects  are  computed  in  the 
frequency  domain  with  the  use  of  the  method  of  moments  (MOM). 
Although  this  method  computes  the  scattered  fields  accurately,  the 
computational  complexity  is  shown1  to  be  0(n3),  where  n  is  the  number 
of  elements  for  which  currents  need  to  be  computed.  Because  of  the 
exponential  nature  of  the  computational  complexity,  as  n  increases,  the 
use  of  MOM  for  computing  scattered  fields  becomes  less  and  less  attrac¬ 
tive.  Other  techniques  for  computing  frequency  domain  scattered  fields 
exist,  namely,  the  fast  multipole  method  (FMM),  multilevel  FMM,  etc.  The 
multilevel  FMM  techniques  have  the  complexity  0(n  log  n).1  Within  the 
realm  of  MOM,  other  techniques  exist  that  take  advantage  of  target 
geometry,  such  as  exploiting  body  of  revolution  (BOR)  symmetry  to 
reduce  the  size  of  matrixes,  which  must  be  inverted.  In  spite  of  great 
advances  in  computational  electromagnetics,  computing  the  scattered 
fields  from  an  arbitrary  target  can  be  prohibitive,  especially  as  the  target 
is  electrically  large.  However,  at  higher  frequencies  for  targets  with 
dimensions  greater  than  2A,  where  A  denotes  the  wavelength,  computing 
scattered  fields  by  physical  optics  (PO)  approximations  yield  relatively 
accurate  results.  Some  of  the  commercially  available  high-frequency 
modeling  packages,  such  as  XPATCH,  model  targets  in  free  space  with  or 
without  conducting  ground  planes.  The  PO  code  that  was  developed 
computes  the  scattered  fields  in  a  lossy  and  dispersive  half-space  via  the 
layered  medium  dyadic  Green's  function  specialized  for  a  half-space.  This 
capability  is  of  interest  for  simulating  targets  in  real  soil.  The  PO  code 
allows  one  to  compute  the  radar  cross  section  (RCS)  and  scattered  fields 
of  an  arbitrary  target  that  is  underground,  above  ground,  or  partially 
buried. 

Section  2  presents  the  basic  formulation  of  the  PO  approximations  and 
the  Green's  functions  used.  Section  3  presents  some  of  the  RCS  character¬ 
istics  of  various  targets,  namely,  a  sphere,  a  vertical  plate,  a  mine,  and  an 
unexploded  ordnance  (UXO),  that  were  computed  with  both  MOM  and 
PO  for  comparison.  Section  3  also  presents  some  of  the  actual  synthetic 
aperture  radar  (SAR)  images  of  UXOs  along  with  the  images  generated 
by  the  simulated  data  by  PO.  These  images  show  that  most  of  the  salient 
features  of  the  targets  are  captured,  which  in  turn,  help  develop  auto¬ 
matic  target  detection  (ATD)  and  automatic  target  recognition  (ATR) 
algorithms. 


1N.  Geng,  A.  Sullivan,  and  L.  Carin,  "Fast  multipole  method  for  scattering  from  3-D  PEC  targets  situated 
in  a  half-space  environment,"  Microwave  Opt.  Tech.  Let.  21,  No.  6  (June  1999),  399-405. 
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2.  Scattering  of  Radiation  by  an  Object  in  Free  Space 


One  assumes  that  the  observations  are  made  in  the  far  field;  that  is, 

(fir  »  1),  where  P  =  cofiue  and  fi  and  e  represent  the  permeability  and 
permittivity  of  free  space,  respectively;  ft  is  the  wave  number;  and  r 
denotes  the  distance  from  the  target  to  the  observation  point.  Figure  1 
shows  the  coordinate  system  used  for  computing  scattered  fields  in  the 
far  field.  I  approximate  the  distance  R  from  any  point  on  the  object  to  the 
observation  point  by 


(r-r'  cosy/  for  phase  variations  , 

~  \r  for  amplitude  variations  , 

where  y/is  the  angle  between  the  vectors  r  and  r'  (vectors  are  denoted  by 
bold  letters).  Then  the  electric  and  magnetic  fields  radiated  by  the  object 
in  the  far  field  are  given  by 


Ea  =  -jco A  (6  and  (j)  components  only)  , 
Hf  =  -jcoF  (6  and  <p  components  only)  , 


where 


A 


lie-iP 

Akt 


N  ,  and 


“■/I  Jsexp (jfir'  cos  y/)  ds  , 


(3) 

(4) 


Figure  1.  Coordinate 
system  for  radiation 
by  an  object. 


R 


where  Js  and  M,  denote  the  electric  and  magnetic  current  densities  in¬ 
duced  on  the  object.  Clearly,  scattered  field  E  can  be  computed  if  Js  is 
known.  To  compute  Js  in  an  object,  one  assumes  a  known  plane  wave  is 
impinging  on  the  object.  PO  approximation  gives  a  relation  between  the 
induced  magnetic  field  and  the  current  density  Jg.  The  following  subsec¬ 
tion  presents  the  relevant  concepts. 


2.1 


Physical  Optics  Equivalent 


For  the  sake  of  completeness,  I  present  in  this  subsection  a  brief  descrip¬ 
tion  of  PO  approximation  theory.  Detailed  explanation  and  derivation  of 
relevant  formulation  can  be  found  in  literature  by  Balanis.2 


Let  us  consider  an  unbounded  medium  with  constituting  parameters  £1 
and  fJv  Let  us  also  assume  that  a  current  source  ]l  radiating  Ex  and  a 
magnetic  current  source  radiating  Hx  are  everywhere  in  the  medium. 
Moreover,  consider  an  imaginary  region  Vy  enclosed  by  a  surface  S1  in 
this  medium  as  shown  in  figure  2.  Assume  that  the  surface  S1  is  replaced 
by  another  medium  with  parameters  e2  and  /L-  Let  us  also  assume  that 
the  sources  and  are  allowed  to  radiate  in  the  presence  of  the  new 
medium  in  Vv  The  total  field  outside  the  region  Vv  produced  by  J1  and 
M:,  is  E  and  H  and  inside  V1  is  Ef  and  Hf. 

The  total  field  outside  the  region  V1  is  equal  to  the  original  field  (E1  and 
Hj)  in  the  absence  of  the  obstacle  plus  a  perturbation  (scattered)  field  (Es, 
Hs),  introduced  by  the  obstacle.  Hence, 


E  =  E:  +  Es  , 

H  =  H1  +  HS  . 


(5) 


In  the  above  equation,  one  can  compute  the  original  fields  E1  and  H1  by 
solving  Maxwell's  equations.  To  compute  Es  and  Hs,  I  introduce  bound¬ 
ary  currents  and  Mp  on  the  surface. 

These  currents  can  be  expressed  in  terms  of  the  difference  in  fields  (Es,  Hs) 
that  are  outside  the  region  Vl  and  the  fields  (Ef,  Hf)  that  are  inside  the 
region  Vy  as  follows: 

Jp  -  ft  x  (Hs  -  Hf)  , 

Mp  =  -«x(Es-Ef)  , 


Figure  2.  Field 
geometry  for 
derivation  of  currents 
on  a  surface  of  an 
object. 
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2C.  A.  Balanis,  Advanced  Engineering  Electromagnetics,  John  Wiley  and  Sons,  Inc.,  New  York,  NY  (1989). 
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where  "x"  denotes  the  cross  product.  Since  the  tangential  components  of 
E  and  H  must  be  continuous  across  the  boundary,  I  get 

Elltan+ESltan=  E*  Itan  x  (E1  +  ES)  =  "  x  E'  , 


H1  Itan  +  RS  Itan  =  ^  I  tan  =>  «  x  (H1  +  Hs)  -  n  X  Hf  , 


(7) 


or 


ES||m-E'ltan=-EllBn=>fiX(ES--E')=-fixEl  - 

Hsl,a„-Ht|Hn=-H1|Mn^flx(Hs-H,)=-fixH1  . 
Substituting  equations  (8)  in  (6),  I  get 


(8) 


Jp  =  -nXH1 
Mp  =  n  xEj  . 


(9) 


Let  us  now  assume  that  the  obstacle  occupying  the  region  is  a  perfect 
electric  conductor  (PEC)  with  conductivity  <7=°°.  Then  the  fields  in  the 
region  Ef  and  Hf  equal  zero  as  shown  in  figure  3.  Over  the  boundary  S  j  of 
the  conductor,  the  total  tangential  components  of  the  E  field  are  equal  to 
zero,  and  the  total  tangential  components  of  the  H  field  are  equal  to  the 
induced  current  density  Jp.  Hence, 

Mp  =  -n  x  (E  -Ef)  =  -it  xE  =  -n  x  (Ej  +  Es)  =  0  , 


or 


-fixEj^nxE5  ,  and  (10) 

Jp=nx(H-Ht)  =  nxH=nx(H1  +  Hs)  , 

or  for  an  infinitely  large  object,  it  can  be  shown  by  image  theory2  that 
H1 1  tan  =  Hs  I  tan.  Hence,  the  PO  for  equation  (10)  becomes 

Jp  =  2n  x  H  j  in  the  lit  region  and 

(11) 

Jp  =  0  in  the  shadow  region  . 


Figure  3.  Field 
geometry  for  physical 
optics  approximation. 
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2C.  A.  Balanis,  Advanced  Engineering  Electromagnetics,  John  Wiley  and  Sons,  Inc.,  New  York,  NY  (1989). 
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Equation  (11)  allows  us  to  compute  the  current  density  Jp  on  the  target 
because  of  Hj,  which  one  assumes  can  be  computed.  Once  we  have 
determined  Jp,  we  can  use  equations  (2)  to  (4)  to  compute  the  field  radi¬ 
ated  by  the  target  with  current  density  Jp.  In  particular,  the  scattered  field 
can  be  computed  once  N  in  equation  (4)  is  determined;  namely. 


N  =  JJ^  J  exp [jpr'  cos  \|/)  ds  =  Jj^  J  exp(//?r  •  r')  ds  ,  (12) 


where  denotes  the  dot  product.  Note  that  J  in  equation  (12)  is  com¬ 
puted  by  equation  (11).  Expanding  equation  (12)  in  terms  of  its  x,  \j,  and  2 
components  of  the  current  J  produces 


N  = 


+  dyJv  +  dzJz)exp(j/3r-r/)ds  . 


(13) 


Using  the  rectangular-to-spherical  component  transformation,  where 


u  ^ 

sin#  cos 0  cos#  cos <p  -sin0 

dr 

a 

— 

sin#  sin  0  cos#  sin  0  cos0 

dg 

V 

S-z 

cos#  -sin#  0 

dp 

we  can  write  the  0  (vertical  polarization)  and  (j>  (horizontal  polarization) 
components  of  N  as 

Ne  =  jj  ( }x  cos  *  cos  0  +  Jy  cos  0  sin  0  -  Jz  sin  0)exp( -r^ds  ,  (15) 

N, p  =  JJ  [Jx  sin  (f)  +  Jy  cos  0)exp(-;/3r  ■  r')ds  .  (16) 

Using  equations  (15)  and  (16)  in  equations  (2)  and  (3)  produces  expres¬ 
sions  for  scattered  field  in  vertical  and  horizontal  polarizations: 

vertical  polarization 

E  q  =  -  e~  ^  ( Jx  cos  ^  cos  <P  +  Jy  cos  ^  sin  <I>-Jz  s'n  0) 

(17) 

x  exp [jP{x' sin  0cos  <p+  y'sin  6>sin  (j)  +  z' cos  0)]  ds 
horizontal  polarization 

(~JX  sin  0  +  /vcos^j  xexp[;/3(x/sin6cos  0  +  y '  sin  0  sin  </)+  z'  cos  6)]ds  (18) 

Note  that  equations  (17)  and  (18)  correspond  to  the  scattered  fields  in  free 
space.  The  following  section  presents  the  half-space  approximations.  We 
can  obtain  these  approximations  by  incorporating  the  half-space  dyadic 
Green's  functions.* 


*T.  Dogaru,  private  notes.  Electrical  Engineering  Dept.,  Duke  University,  Durham,  NC. 
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2.2  Half-Space  Approximations 

Consider  two  layers  with  an  xy-plane  with  z  =  0  as  the  boundary  between 
the  two  as  shown  in  figure  4.  The  medium  in  layer  1  has  the  permittivity 
e1  and  permeability  nx  while  the  medium  in  layer  2  has  e2  and  y2.  Often 
the  medium  in  layer  1  is  air;  thus,  nx  =  jiQ  =  47rlCT7  H/m  and  £1  =  £  0  = 
10'9/367r  F/m.  Define  the  following  functions: 


Pi  =  > 


Pi  =  (Offhei  , 

expx  =  exp[//?1(x'sin0cos0  + y'sinflsin^  +  z'cos#)]  , 
exp  2  =  exp[;  fi  1  (x '  sin  6  cos  <p  +  y '  sin  0  sin  (f>  -  z '  cos  0)]  , 

exp3  =  exp[/(/?1(x'sin  0cos0  +  y/sin0sin^)  +  z'\f P\- (/i}Sin0)2 )]  , 


T  =T  = 

XX  1yy 


IfilCOSd 


T  = 

1  22 


Pxcos6  +  \f p*  -  (/?! sin0)2 
2/?x  cos  0 


T  = 

1zx 


j^cos0  +  V  ^-(ftsin©)2 
^cosg-V^-(^sing)7  ^cosg_^|T 


A 


(/^sinflV 


^cos«  +  >/^-(/3i8in(?7'  Pl^se +MIiP^? 


cos  0  cos  (j) 


sin  6 


sin  0 

1 2y  -  J  cos  (/>  ’ 

»  jSiC  os0-y//^-(/?1sin0)2 

xx  w~  plCose+<M -  iPiSine)2  ' 

r,  _D  )8iCOsg-//3i-(/?iSin0)T 

a~  hcose  +  Jpl-faime)2  ' 


R  =T 
iXzx  1  zx 


R zy  ~  R-zx 


sin  0 
cos  (f) 
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Figure  4.  Geometry 
used  for  layered 
medium. 
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The  following  are  the  expressions  for  Ee  and  E.  for  the  target  in  layer  1: 


Ee  =  ~^KrLe~^l'  jis  ^C°S  6C°S  ^eXPl  +  Rxx exp2)  “ sin  9Rzx exP 2] 

+  /y[cos  6  sin  <p(exp1  +  Rxyxj  exp2)  -  sin  6RZX/  exp2]  -  Jz  sin  0(expj  +  Rzz  exp2))  ds  , 

Et>  =  ~  /f^e_//?irjJs  (~Ix  sin^eXPl  +  Rxx  exp2)  +  h  cos^(exPl  +  Ryy  exP 2)) ds  • 

Expressions  for  Ee  and  for  the  target  in  layer  2  are  given  by 

Ee  =  ~  ^~e~iPv  j\s  (/x(cos  9  cos  <t>Txx  -  sin  6  Tzx)ex p3  ds 

+  /v(cos  0  sin  0  Tyy  -  sin  6  T2J/)exp3  -  Jz  sin  6  Tzz  exp3)  ds  , 


(19) 

(20) 

(21) 


-iPv  JIh  x  sin  0  Txx  exp  3  +  Jy  cos  0  TXJXJ  exp3)  ds  . 


Equations  (19)  to  (22)  can  be  used  for  computing  the  scattered  fields  by  a 
target. 


2.3  Implementation  of  PO  Approximations  in  Half-Space 
for  Computation  of  Scattered  Fields 

To  compute  the  scattered  fields  near  a  target,  one  must  first  partition  the 
target  into  small  triangles  with  each  side  of  the  triangle  approximately 
0.1A,  where  X  denotes  the  wavelength.  For  each  triangle  patch,  an  out¬ 
ward  normal  and  its  centroid  are  computed.  This  normal  is  used  to 
compute  the  current  induced  by  a  magnetic  field  in  the  patch.  The  cen¬ 
troid  becomes  the  reference  point  in  the  patch  for  computing  electric  and 
magnetic  fields.  The  magnetic  field  in  a  patch  is  then  used  to  determine 
the  current  with  the  use  of  PO  approximation.  Once  the  current  in  a  patch 
is  computed,  it  is  then  used  to  compute  the  electrical  field  radiated  (scat¬ 
tered  field)  by  the  patch.  The  sum  total  of  all  the  fields  radiated  by  all  the 
patches  constitutes  the  scattered  field  by  the  target. 


Now  to  compute  the  induced  current  caused  by  a  magnetic  field,  let  us 
assume  that  an  incident  plane  wave  where  the  electric  field  has  a  magni¬ 
tude  E0  =  1  is  impinging  on  the  target.  Let  the  direction  of  the  incident 
plane  wave  be  given  by  9  and  (f>;  then  the  unit-transmitted  vector  in 
rectangular-coordinate  system  is  given  by 

layer  1:  tx  =  -  [sin  9  cos  (j>,  sin  9  sin  0,  cos  9 ] 
r  =  k|tx, 

and  the  reflected  vector  is  rj  =  -[sin  9  cos  (j),  sin  6  sin  (f>,  -cos  6],  and 
layer  2:  tx  =  -  [sin  6t  cos  <j),  sin  6t  sin  0,  cos  9t], 

where  6t  denotes  the  angle  of  refraction  in  medium  2  and  is  determined 
by  Snell's  law:  sin  6t  =  jS1  sin  6/  p2  ■ 

Figure  5  shows  the  transmitted,  reflected,  and  refracted  rays  (vectors)  and 
the  corresponding  angles. 

Field  computations  in  layer  1: 

The  electric  fields  in  vertical  and  horizontal  polarizations  impinging  on  a 
patch  are  given  by 

transmitted  field 

Ey  =  Eo[cos0cos0,  cos 0  sin 0,  -sin0]exp[-/yS1(tx  •  P,)]  , 

(23) 

Eh  =  E0[-  sin  (j),  cos  (p,  0]exp\-i/d1(tx  ■  P/)j  , 

where  P ;  =  [x-,  y  -,  z-]  denotes  the  coordinates  of  the  centroid  of  an  zth 
triangular  patch  of  the  target. 


Figure  5.  Ray 
dynamics. 
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reflected  field 

E y  =  (Eort,)[cos0  cos (f>,  cos 0  sin 0,  sin0]exp[-//?1(ri  -  P;)]  , 

eH  =  Ml~  sin  <t>,  cos  (j),  0]exp[-;^1(r1  •  P;)]  , 

where  the  reflection  coefficients  rv  and  rh  are  given  by 

77  2  cos 0t-  r/1  cos 0 
Ty  ~  i]2  cos6t  +  cosO  ' 


(24) 


and 


?72cos0-  rj^  cos 6t 
r,z  -  ri2  cosd  +  77 2  cos  0^ 


Vl  = 


772  = 


0\ 


intrinsic  impedance  of  medium  1  , 
intrinsic  impedance  of  medium  2  . 


Then  the  magnetic  fields  can  be  expressed  in  terms  of  the  electric  fields  as 

Hy  =  7^'(txxEv)  ' 


hh  -  x  eh)  / 
Hrv  =  ±(rlxVv)  , 


(25) 


where  a  x  b  represents  the  cross  product  between  two  vectors  a  and  b. 

Computation  of  currents  in  a  triangular  patch: 

Having  determined  the  magnetic  fields  in  a  triangular  patch,  I  can  now 
compute  the  currents  using  the  PO  approximation  given  by  equation  (11) 
as  follows: 


Jy-2(fipXHy) 


Jil  -  2(fipX  Hpjr) 


Jy  =  2(fip  X  Hy) 


(26) 


Jpf-2(npxHy  , 

where  np  denotes  the  outward  normal  to  the  triangular  patch  P.  Note 
that  the  normal  rip  and  the  Hy  should  correspond  to  the  same  triangular 
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patch.  Now  I  can  compute  the  scattered  fields  using  equations  (19)  and 
(20).  Note  that  equation  (19)  gives  the  scattered  field  in  vertical  polariza¬ 
tion.  By  using  Jy,  the  current  caused  by  the  vertically  polarized  transmit¬ 
ted  field  in  equations  (19)  and  (20),  I  can  obtain  scattered  fields  in  vertical 
transmit,  vertical  receive  (VV)  and  vertical  transmit,  horizontal  receive 
(VH)  polarizations.  Similarly,  using  in  equations  (19)  and  (20),  I  can 
obtain  scattered  fields  in  HV  and  HH  polarizations: 

Eyy  =  -  Uvx[cos  9  cos  0(exPl  +  Rxx  exp2)  -  sin  9RZX  exp2] 

(27) 

+  Jvy[c os  0  sin  ^(exp1  +  Ryy  exp2)  -  sin  9Rzy  exp2]  -  }Vz  sin  d(exp]  +  Rzz  exp 2))ds  , 

Evh  =  - Hv*sin0(exPi  +  Rxx exP2)  +  7vyCos0(exPl  +  Ryy exp2) ) ds  ,  (28) 


E 


=  “  ^HxtC°S  6  COS  ^6Xp  1  +  R*x  exP2)  -  sin  dRzx  exPl] 

+  Jhi/[c os  0  sin  0(exp2  +  Ryy  exp2j  -  sin  9RZy  exp2j-/Wzsin^(exp1  +  Rzz  exp2)  ds  , 


(29) 


aHH  =  - 


mh 

4  7tR  1 


Wjls  Hhx  sin  0(exPi  +  RxxexPl)  +  Jhij  cos  0|  exp  \  +  Ryy  exp2) )  ds  ,  (30) 


where  Jv  =  [JVx,  JVy,  JVz\  and  R  is  the  distance  between  the  observation 
point  and  the  target.  Similarly,  E  yy,  E  yW,  E  hv>  and  ^HH  can  be  com¬ 
puted  with  similar  equations  (eq.  (27)  to  (30)),  where  Jy  and  are  re¬ 
placed  by  J  v  and  J  fo,  respectively.  These  scattered  electric  fields  are 
computed  for  only  those  triangle  patches  illuminated  by  either  tx  or  q. 
The  sum  of  Eyy  +  Eyy  for  all  the  illuminated  patches  gives  the  scattered 
electric  field  in  W  polarization.  Similarly,  other  summations  give  the 
field  in  VH,  HV,  and  HH  polarizations. 

Computations  in  layer  2: 

Field  computations  in  layer  2  are  as  follows: 


Ey  =  (E0  ■  Ty)[cos0f  cos 0,  cos  9t  sin  0,-sin  9t]exp[-j(31(lx  ■  Pz)]  , 

Eh  =  (E0th)[- sin^' cos^  °]exP  [-jP l(tx  p/)]  / 


(31) 


where 


tx  =  -[sin  9t  cos  0,  sin  9t  sin  0,  cos  9t] , 
and  the  transmission  coefficients  Ty  and  TH  are  given  by 
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T  2?72  cos  6 

V  ?12COS0f+  7]  I  cos  6 

X  2772  cos  g 

H  772  cos0  +  rh  cos  0t 

Notice  that  no  reflected  fields  are  in  layer  2.  As  in  layer  1,  magnetic  fields 
and  currents  in  each  triangular  patch  are  computed,  and  the  scattered 
fields  are  computed  with  equations  (27)  to  (30)  and  appropriate  currents. 
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3.  Validation  of  PO  Code 


A  PO  code  was  implemented  in  a  matrix  laboratory  (MATLAB)  environ¬ 
ment.  To  verify  the  accuracy  of  the  code,  I  created  several  targets  and 
computed  their  RCSs  using  both  PO  code  and  MOM  code.  The  MOM 
code  is  developed  by  Duke  University.  The  comparison  results  are  pre¬ 
sented  below. 

Vertical  plate:  A  1-m2  plate  was  created  as  shown  in  figure  6  and  placed 
0.15  m  above  the  Eglin  soil  with  permittivity  of  £,.  =  5.0  -jO.O  and  conduc¬ 
tivity  of  a  =  0.003.  For  the  computation  of  the  scattered  fields  at  various 
frequencies,  an  incident  field  at  d  =  45°  and  0  =  0°  was  used.  The  plate 
was  oriented  to  get  the  dihedral  effect  from  the  ground.  I  computed  the 
RCS  characteristics  of  the  plate  from  25  to  2000  MHz.  The  results  are 
shown  in  figures  7  and  8  for  VV  and  HH  polarizations,  respectively. 

Figures  7  and  8  show  that  PO  output  is  fairly  accurate  compared  to  that  of 
the  actual  output  computed  by  MOM  at  frequencies  600  MHz  and  higher. 
In  other  words,  when  the  target  was  greater  than  or  equal  to  2A,  I  ob¬ 
tained  good  correlation  with  the  MOM  simulations.  Why  the  PO  output 
for  HH  polarization  had  better  correlation  at  even  lower  frequencies  with 
MOM  output  compared  to  W  polarization  is  not  clear. 

The  next  object  I  used  for  simulations  and  comparison  was  a  horizontal 
cylinder  with  a  spherical  head  on  one  side  as  shown  in  figure  9.  The 
dimensions  of  the  cylinder  were  length  =  0.58  m  and  diameter  =  0.078  m. 

I  performed  two  simulations,  namely,  a  cylinder  buried  1  in.  below  Yuma 
soil  and  a  cylinder  placed  1  in.  above  Eglin  soil.  The  results  of  the  simula¬ 
tions  are  shown  in  figures  10  to  15  for  a  target  buried  1  in.  below  the 
surface  in  Yuma  soil.  PO  output  is  plotted  by  a  solid  line  and  MOM 
output  is  plotted  by  a  dotted  line.  Notice  from  these  figures  that  the  RCS 
computed  for  the  cylinder  using  the  PO  code  is  accurate.  The  largest 
difference  between  the  MOM  and  PO  code  is  observed  in  figure  15,  where 
<t>  =  -90°;  at  that  angle,  it  is  known2  that  PO  does  not  do  well.  Figures  16  to 
21  show  the  RCS  characteristics  of  a  cylinder  placed  1  in.  above  the 
surface  of  Eglin  soil. 

First,  I  simulated  the  scattering  fields  for  the  cylinder  shown  in  figure  9 
for  all  azimuthal  angles  (p  =  0°  to  360°  and  generated  SAR  images  using 
the  back-projection  algorithm.  The  SAR  images  generated  using  MOM 
and  PO  simulated  data  are  shown  in  figure  22  for  a  cylinder  buried  1  in. 
below  the  surface  of  Yuma  soil.  From  the  SAR  images,  I  found  that  an 
excellent  correlation  exists  between  the  MOM  and  PO  simulated  SAR 
images  in  terms  of  resolution  and  that  both  the  images  have  similar 
magnitudes.  Next,  I  used  a  UXO  (shown  in  fig.  23),  which  was  buried  in 
Yuma  soil.  I  made  actual  measurements  of  the  UXO  using  a  radar.  Its  SAR 
images  and  the  PO  simulated  images  for  the  same  UXO  are  presented  in 
figure  24.  It  shows  that  the  SAR  images  using  PO  simulated  data  can 
capture  all  the  features  of  the  UXO.  This  helps  in  developing  ATD  and 
ATR  algorithms  for  various  targets. 

2C.  A.  Balanis,  Advanced  Engineering  Electromagnetics,  John  Wiley  and  Sons,  Inc.,  New  York,  NY  (1989). 
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Figure  14.  RCS  VV  of 
a  cylinder,  0  =  60°, 

0  =  _90°  (axial). 


Figure  15.  RCS  HH  of 
a  cylinder,  9  =  60°, 

0  =  - 90°  (axial). 
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Figure  16.  RCS  VV  of 
a  cylinder  above 
Eglin  soil,  6  =  60°, 

0  =  0°  (broadside). 


Figure  17.  RCS  HH  of 
a  cylinder  above 
Eglin  soil,  6  =  60°, 
0=0°  (broadside). 
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Figure  18.  RCS  VV  of 
a  cylinder  above 
Eglin  soil,  0  =  60°, 

(p=  45°  (off  axis). 


Figure  19.  RCS  HH  of 
a  cylinder  above 
Eglin  soil,  9  =  60°, 

(J)  =  45°  (off  axis). 
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4.  Remarks 


One  of  the  drawbacks  in  the  current  PO  code  is  that  normals  of  various 
patches  are  being  used  to  determine  whether  a  patch  is  illuminated  by  the 
transmitted /reflected  ray.  However,  it  does  not  determine  whether  a 
patch  is  blocked  by  some  part  of  the  target.  For  some  simple  targets,  this 
is  not  a  problem.  If  the  target  is  complex  and  has  many  corners  and 
projections,  the  PO  code  may  not  give  accurate  results  unless  the  above- 
mentioned  problem  is  solved. 
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